Phonon Linewidths and Electron Phonon Coupling in Nanotubes 
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We prove that Electron-phonon coupling (EPC) is the major source of broadening for the Raman 
G and G~ peaks in graphite and metallic nanotubes. This allows us to directly measure the optical- 
phonon EPCs from the G and G~ linewidths. The experimental EPCs compare extremely well with 
those from density functional theory. We show that the EPC explains the difference in the Raman 
spectra of metallic and semiconducting nanotubes and their dependence on tube diameter. We 
dismiss the common assignment of the G~ peak in metallic nanotubes to a Fano resonance between 
phonons and plasmons. We assign the G + and G~ peaks to TO (tangential) and LO (axial) modes. 

PACS numbers: 73.63.Fg, 63.20.Kr, 73.22.-f 
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Electron-phonon coupling (EPC) is a key physical pa- 
rameter in nanotubes. Ballistic transport, supercon- 
ductivity, excited state dynamics, Raman spectra and 
phonon dispersions all fundamentally depend on it. In 
particular, the optical phonons EPC are extremely rele- 
vant since electron scattering by optical phonons sets the 
ultimate limit to high field ballistic transport Q, 0, 0, 
0, . Furthermore, they play a key role in defining the 
phonon dispersions and the Raman spectra of metal- 
lic and semiconducting single wall nanotubes (SWNT), 
which show different features [1 B II ES El El El 
El EH EU El E3 • Several theoretical and experimen- 
tal investigations of acoustic phonons EPC have been 
published (see, e.g., Ref. However, only tight- 

binding calculations of optical phonons EPC were per- 
formed, with contrasting results [1 E M Ell El 12^. 
More crucially, no direct measurement of optical phonons 
EPC has been reported to validate these calculations. 

In this Letter we prove that the optical phonons EPC 
are the major source of broadening for the Raman G and 
G~ peaks in graphite and metallic SWNTs. We show 
that the experimental Raman linewidths provide a direct 
EPC measurement. The EPC are also responsible for the 
the G+ and G" splitting in metallic SWNTs. This allows 
us to unambiguously assign the G + and G~ peaks to TO 
(tangential) and LO (axial ) m odes, in contrast to what 
often done l^lil llHllulTilTll^ . We dismiss the common 
assignment of the G~ peak in metallic SWNTs to a Fano 
resonance between phonons and plasmons |ToL Ell IT3. Ip3 . 
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In a perfect crystal, the linewidth 7 of a phonon is 
determined by its interaction with other elementary ex- 
citations. Usually, 7 = 7 a "+7 BP , where 7°™ is due to the 
interaction with other phonons and r ) EP with electron- 
hole pairs. 7 a ™ is determined by anharmonic terms in 
the interatomic potential and is always there. r y EP is 
determined by the EPC and is present only in systems 
where the electron gap is zero. If the anharmonic contri- 
bution 7 a " is negligible or otherwise known, measuring 
the linewidth is the simplest way to determine the EPC. 



A phonon is described by a wave vector q, branch in- 
dex -q and frequency w q7? . We consider a mean-field sin- 
gle particle formalism, such as density functional theory 
(DFT) or Hartree-Fock. The EPC contribution to 7 qj) is 
given by the Fermi golden rule |2(| : 
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^ l5(k+q)j,ki| 2 [/ki - /(k+q)j] 
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where the sum is on the electron vectors k and bands i 
and j, iVk is the number of k vectors, /ki is the occupa- 
tion of the electron state |k,z), with energy e^i, S is the 
Dirac distribution. <7(k+ q )j,ki = D 



(k+q)j,k 



where M is the atomic mass. -D(k+q)j.ki = \ K ~l~ 
q, j|AVq^|k, i) and AVq^ is the potential derivative with 
respect to the phonon displacement. D is the EPC. 

The electron states contributing to the sum in Eq. 
are selected by the energy conservation condition eti + 
hjj^ — e(k+q)j- Also, the state ki has to be occupied 
and (k + q)j empty, so that the term [/kj — /(k+q)j] 7^ 0. 
Thus, only electrons in the vicinity of the Fermi level 
contribute to r y EP . In insulating and semiconducting sys- 
tems 7 BP = 0. In general, a precise estimate of r y EP from 
Eq. n is possible only after an accurate determination of 
the Fermi surface. However, graphite and SWNTs are 
very fortunate cases. Thanks to their particular band 
structure, j EP is given by a simple analytic formula. 

In general, the EPC determines both phonon dis- 
persions and linewidths. We first consider the case of 
graphite and show that both dispersion and linewidths 
give a direct measure of the EPC at T. 

The electron bands of graphite are well described by 
those of a two dimensional graphene sheet. In graphene, 
the gap is zero for the tt bands at the two equivalent K 
and K' = 2K points of the Brillouin zone. We define 
{D^)f = Y^l j I^Kj,Kj | 2 /4, where the sum is on the two 
degenerate it bands at the Fermi level €f- We consider 
the EPC relative to the E 2g phonon at T. The r-E 2s 
mode is doubly degenerate and consists of an antiphase 
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FIG. 1: Graphite phonon dispersion of the highest optical 
branch near T. Dots are inelastic X-ray measurements from 
Ref. |2^|. The line is a quadratic fit. 



in-plane motion. For a small non zero q near T, this 
splits into a quasi longitudinal (LO) and quasi transverse 
(TO) branch, corresponding to an atomic motion parallel 
and perpendicular to q. From DFT calculations we get 
(D 2 ) F =45.60(eV/A) 2 for both LO and TO modes @. 

We have previously shown that graphite has two Kohn 
anomalies in the phonon dispersions for the T-^g and K- 
modes 6] . Due to the anomaly, the dispersion near F 
of the E2 S -LO mode is almost linear, with slope S 1 - |(|: 



sh° 



V3hal 
8Mu T /3 



(D 2 t )f, 



(2) 



where ao=2.46 A is the graphite lattice spacing, 
(3 =5.52 A eV is the slope of the electron bands near 
ep, M is the carbon atomic mass and wp is the fre- 
quency of the E 2s phonon (huj r =196.0 meV). Eq.|2shows 
that (-Dp) f can be directly measured from the exper- 
imental Sp . The phonons around V have been mea- 
sured by several groups with close agreement. From a 
quadratic fit to the most recent data of ref. [2?] (Fig. 
we get Sh° = 133 cm^'A. From Eq.Qwe have (D^) F = 
39 (eV/A) 2 , in good agreement with DFT. 

We now show that the Full Width at Half Maximum 
(FWHM) of the graphite G peak gives another indepen- 
dent EPC measurement. The G peak of graphite is due 
to the r-E 2g phonon [28| . We use Eq. to compute the 
width, 7p P , for this mode. Close to K, we assume the 
7r bands dispersion to be conic from the Fermi level ep, 
with slope (3 (Fig.EJi) H^. For both LO and TO modes: 



V3a 2 h 2 
AM/3 2 
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We then measure FWHM(G) for a single-crystal graphite 
similar to that of Ref. [2jj . Its Raman spectrum does not 
show any D peak (Reich in and E3l), thus we exclude 
extra broadening due to disorder [3l|. By fitting the 
G peak with a Lorentzian we get FWHM(G) = 13 cm -1 . 
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FIG. 2: Electron bands around K in graphene (a) and in a 
metallic tube (b). Shaded area is the graphene Brillouin zone. 
Dashed arrows: decay processes for a T phonon. 



Temperature dependent measurements show no increase 
of FWHM(G) in the 2K-900K range EI Our Raman 
spectrometer resolution is ~ 1.5cm -1 [30|. We thus as- 
sume an anharmonic contribution lower than the spec- 
tral resolution Jq 1 < 1.5cm _1 . Thus, we estimate 7p P ~ 
11.5 cm" 1 . Then, from Eq.El{£) r ) F - 47 (eV/A) 2 . This 
compares very well with DFT, again confirming that Jq 1 
is small. 

Finally, near T the conservation of the energy and mo- 
mentum in Eq. ^ implies: 



q > hior/ (3. 



(4) 



This condition is satisfied by the E2 9 phonon, involved in 
the Raman G peak of graphite. On the other hand, the 
double resonant mode close to T, which gives the D' peak 
at~1615 cm -1 , does not satisfy it. D' is indeed sharper 
than the G peak . 

In summary, we presented two independent measure- 
ments of the graphite r-E 2£( EPC. These are consistent 
with each other and with DFT. We now consider SWNTs. 

The G peak of SWNTs can be fit with only two com- 
ponents, G + and G~ 0. Semiconducting SWNTs have 
sharp G + and G~ whilst metallic SWNTs have a broad 

downshifted g- 0li,liEElElIilElElElE3E3 

The G~ band shows a strong diameter dependence, be- 
ing lower in frequency for smaller diameters |l l| . This 
suggested its attribution to a tangential mode, whose 
circumferential atoms displacements would be most af- 
fected by a variation in diameter 0. Thus, the G+ 
and G - peaks are often assigned to LO (axial) and TO 
(tangential) modes 0, H El H El 

Conflicting reports exist on the presence and rela- 
tive intensity of the G~ band in isolated versus bundled 
metallic tubes. It has been claimed that this peak is as 
intense in isolated SWNTs as in bundles |llHl7l ll8 l : that 
it is smaller [r2lll4lll^ |: or that it can be absent [24(. The 
G~ peak is also thought to represent a Fano resonance 
due to phonon coupling with plasmons [ToL fl^. IT3L 12^ . l25|| . 
Such phonon-plasmon coupling would either need |lj| 
or not need jlClL a finite phonon wavevector for its 
activation. The theory of Refs. predicts the 

phonon-plasmon peak to be intrinsic in single SWNT, 
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FIG. 3: (color online) Black points: experimental G + and 
G" in metallic SWNTs vs. diameter. Red points: DFT for 
armchair [from (3,3) to (18,18)], zigzag [from (9,0) to (30,0)] 
and chiral tubes [ (5,2), (12,3), (16,1), (16,10), (20,14)]. Lines: 
fit of Eq. 5 to the experimental and DFT data. 
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FIG. 4: Experimental FWHM(G") in metallic SWNTs vs 
ameter. Solid line: DFT. Dotted line: fit of FWHM 
from Eq.Eland -y an = 10cm -1 . 
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in contrast with [24J. On the other hand, the theory in 
Ref. E3 requires several tubes (>20) in a bundle in or- 
der to observe a significant G~ intensity, in contrast with 
the experimental observation that bundles with very few 
metallic tubes show a significant G~ El El El Hi 
Ref. E3 a l so predicts a G~ upshift with number of 
tubes in the bundle, in contrast with Ref. [24], which 
shows a downshift, and with Refs. 0, 0, ^3, which 
show that the G~ position depends on the tube diameter 
and not bundle size. Finally, the G _ position predicted 
by JrIEeI is at least 200 cm -1 lower than that measured 
0llliE3ElElElElElE3E3 Thus, all the pro- 
posed theories for phonon-plasmon coupling 0, 0, |2j| 
are very qualitative, require the guess of several physical 
quantities, and fail to predict in a precise, quantitative, 
parameter-free way the observed line-shapes and their 
dependence on the SWNT diameter. 

We now show that, like in graphite, the EPC per- 
se gives the main contribution to the G _ position and 



FWHM, even in the absence of phonon-plasmon cou- 
pling. Surprisingly, this has not been considered so far. 
However, it is clear that, if phonons do not couple to 
electrons, they certainly cannot couple to plasmons. 

We compute the phonon frequencies for several metal- 
lic SWNTs using DFT (details are in Refs. HH3). For 
all metallic tubes, of any chirality, we get a splitting of 
the modes corresponding to the graphite r-E2 g into quasi 
transverse (TO) and a quasi longitudinal (LO), corre- 
sponding to an atomic motion tangential to the SWNT 
axis or parallel to it. The LO frequency is smaller than 
the TO and has a strong diameter dependence(Fig. |3J. 
This is the opposite of what originally proposed, i.e. that 
a tangential mode should have a much stronger diame- 
ter dependence than an axial one @, H 0, El [3 E3] • 
This counter-intuitive result can be understood only con- 
sidering the presence of a Kohn anomaly in the phonon 
dispersion of metallic SWNTs 0, S Q. The Kohn 
anomaly lowers the frequency of phonons having a sig- 
nificant EPC. While the LO EPC is large, the TO EPC 
is zero. Thus, the LO frequency is smaller than the TO. 
Furthermore, from the dependence of the EPC on the 
tube diameter d we get the analytic result |33| : 



(5) 



with Cm a constant. We thus assign the G~ peak to 
the LO (axial) mode and the G + to the TO (tangential) 
mode. Fig. [21 compares our DFT calculations with ex- 
periments. The agreement is excellent, considering that 
no fitting parameters are used. If we fit Eq. Oto Fig. 01 
we get Ca/=1.52x10 5 cm" 2 nm; ujto— 1591 cm -1 . Note 
that a 1/d 4 dependence of was previously suggested 
based on phonon-plasmon coupling [nj, [lj . 

Now we use Eq. ^ to derive the EPC contribution to 
FWHM(G+) and FWHM(G~) in metallic SWNTs. The 
EPC of a SWNT can be obtained from the graphite EPC 
(Z?p)f via zone-folding (valid for d >0.8 nm, i.e. for the 
vast majority of SWNTs used in experiments and de- 
vices) 3]. Four scattering processes are involved (Fig. [2J- 
b) and the LO and TO linewidths are: 



7r-LO 



(D 2 r ) 



nMujrP d 



7^x0 = 0. 



(6) 



Eq.[|jlis a key result. It shows that the EPC contributes to 
the linewidth only for the LO mode in metallic SWNTs. 
For semiconducting SWNTs the EPC contribution is zero 
for both the TO and LO modes, since the gap does not 
allow to satisfy the energy conservation in Eq. ^ 

This confirms our assignment of the G~ and G + Ra- 
man peaks to the LO and TO modes. Experimentally, in 
semiconducting SWNTs FWHM(G+) and FWHM(G~) 
are similar (both are ~ 10 cm -1 ), while in metallic tubes 
FWHM(G"b: 60 cm -1 and FWHM(G+)~ 10 cm -1 , for 
a given d |ll E3- The large FWHM(G~) in metallic 
SWNTs is due the large j EP and the small FWHM(G+) 
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is entirely anharmonic in origin |35j. Thus, even if 
FWHM(G) in graphite and FWHM(G+) in SWNTs are 
similar, their origin is totally different. This also explains 
why a FWHM(G+)<FWHM(G) can be seen [TJ. 

Eq. explains the 1/d dependence of FWHM(G~) in 
metallic SWNTs Then, using Eq.© (D^>jp can be 
directly fit from the experimental FWHM(G") (Fig.EJ. 
We find (Df~)p = 37(eV/A) 2 , again in agreement with 
DFT. Note that, while the DFT (D%) F is that of planar 
graphene, the fitted {D^)f is obtained from measure- 
ments of SWNTs. Thus, the good agreement between 
the two values is a direct verification that the SWNT 
EPC can be obtained by folding the graphene EPC [3. 

The most common process involved in Raman scatter- 
ing is single resonance. Double resonance is necessary to 
explain the activation of otherwise inactive phonons away 
from r, such as the D peak 0,H3 It has been suggested 
that even the G + and G~ peaks in SWNTs are always 
double resonant 0,|3(|. However, if the laser excitation 
energy satisfies single resonance conditions, the intensity 
of single resonance peaks is expected to be dominant in 
the Raman spectrum |38| . Double resonance can only be 
relevant for higher excitation energies. The condition set 
by Eq. 01 must also hold for SWNTs in order to have a 
significant EPC contribution to the linewidth. Eq.0]can 
only be satisfied by phonons with wavevector too small to 
be double resonant 0,^3- Thus, in double resonance, 
the G _ peak should be much narrower than experimen- 
tally observed. Interestingly, it has been reported that 
the broad G~ peak disappears by increasing the excita- 
tion energy while measuring the Raman spectrum of a 
metallic SWNT, and a sharper one appears Q|. We ex- 
plain this as a transition from single to double resonance. 

Finally, in Ref. [j| we computed the zero-temperature 
life-time r of a conduction electron in a SWNT, due to 
phonon scattering. If phonons are thermalized, t deter- 
mines the electron mean free path Iq = fir/H at zero 
temperature. The phonon linewidths can be expressed 
as a function of r as: 



iT-LO 



LO 



4/rp_ TO , 



(7) 



where bs(fs) indicate back (forward)-scattering [3[. Such 
simple relations might seem surprising. In fact, the two 
quantities 7 -1 and r describe two distinct phenomena: 
7 _1 is the life-time of a phonon (Eq.[T} and t is the life- 
time of an electron (Eq. 1 of Ref. |3j). In general, one 
does not expect a simple relation between r y~ 1 and r. 
However, due to the low dimensionality, in SWNTs r y EP 
determines directly r (Eq. UJ. Thus, the measured r y EP 
are also a direct measurement of Iq for metallic SWNTs. 
Since the measured j EP is in agreement with DFT cal- 
culations, the results of Ref. Q are further confirmed. 

In conclusion, we presented a set of simple formulas 
(Eqs. El EI), with the graphite EPC as the only fit 
parameter. Remarkably, these formulas quantitatively 
explain a string of experiments, ranging from phonon 



slopes and FWHM(G) in graphite to G peak position 
and FWHM(G") diameter dependence in SWNTs. Fit- 
ting a wide set of independent data we obtain the same 
EPCs ±10%. These experimental EPC are in excellent 
agreement with, and validate, the DFT approach. 
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IHP-HPMT-CT-2000-00209 and CANAPE, from EP- 
SRC GR/S97613 and The Royal Society is acknowledged. 
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